Drell-Yan process at forward rapidity at the LHC 



O 

(N 

< 



o 
o 



- 1—1 

X 



Krzysztof Golec-Biernat, 1 ' 2 B Emilia Lewandowska, 2: and Anna M. Stasto 3 ' 4, 2, 

'institute of Physics, University of Rzeszow, Rzeszow, Poland 
^Institute of Nuclear Physics Polish Academy of Sciences, Cracow, Poland 
s Penn State University, Physics Department, University Park, PA 16802, USA 
4 RIKEN Center, Brookhaven National Laboratory, Upton, NY 11973, USA 

We analyze the Drell-Yan lepton pair production at forward rapidity at the Large Hadron Collider. 
Using the dipole framework for the computation of the cross section we find a significant suppression 
in comparison to the collinear factorization formula due to saturation effects in the dipole cross 
section. We develop a twist expansion in powers of Q 2 S /M 2 where Q s is the saturation scale and 
M the invariant mass of the produced lepton pair. For the nominal LHC energy the leading twist 
description is sufficient down to masses of 6 GeV. Below that value the higher twist terms give a 
significant contribution. 



I. INTRODUCTION 



The Large Hadron Collider (LHC) opens a new kinematic regime at high energies. In this regime QCD evolution 
leads to the fast growth of the gluon density. At these high densities it is possible that the novel phenomena related 
to the nonlinear dynamics of the gluon fields will occur. Drell-Yan production is a unique process which offers high 
' sensitivity to the parton distribution in the hadron. It is one of the few processes in hadron-hadron collisions where 
^p" 1 , the collinear factorization has been rigorously proven Within this framework, the NLO calculations have been 

O : performed in [1-0], and later on up to NNLO accuracy in [8l-[l0j. In the collinear factorization approach the DY 
process is viewed as the fusion of the quark and antiquark which produces a virtual (timelike) photon. 

One can view the same process alternatively in the rest frame of one of the hadrons. The quark (typically valence 
one) from the fast incoming projectile interacts with the color field of the target hadron, and emits the virtual photon. 
t-H ■ The photon then decays producing lepton pair which moves into the region of forward rapidity (with respect to the 
incoming projectile). The original description of the DY in the dipole picture has been proposed in [ll|, [l2| with 
details of the calculations presented in [13[ . This process has been also reexamined in [14( and later on in 
within the framework of the Color Glass Condensate. This formulation is applicable to the very forward production, 
when the fractions of the longitudinal momenta of the incoming partons are very different (see discussion in the 
. ■ next section). This approach has been very useful as one can easily incorporate the higher twist effects due to the 
multiple scattering of dipole off the target field. When the energy is high, one of the momentum fractions of the 
incoming partons is very small, and as a result one is probing potentially dense gluon fields in the target. This means 
that the multiple scatterings have to be taken into account. Consequently the parton evolution should be modified 
by inclusion of higher twist terms, i.e. terms which are subleading in the expansion of Q s (x) 2 /M 2 where Q s (x) is 
the saturation scale characterizing the dense gluon field in the target, and M 2 is the invariant mass squared of the 
produced Drell-Yan pair. 

The twist expansion method for analyzing the deep inelastic process has been constructed in [17| for the case of 
the dipole model based on the idea presented in [18(. There, a systematic expansion in powers of Q s (x) 2 /Q 2 was 
performed and terms up to twist 4 were extracted analytically both for the longitudinal and transverse structure 
function. The analysis was further extended to the dipole model which included the DGLAP evolution corrections 

In this paper we develop a twist expansion for the Drell-Yan process in a similar approach to that presented in [171 ] 
for DIS. Unlike the DIS case however, due to the presence of the second hadron (projectile) the expansion involves 
additional dependence on the structure function of that fast moving hadron. As a result the obtained expressions 
still contain the integrals over the longitudinal momentum fractions which involve the structure function. Using this 
formalism we evaluate the twist 2 and twist 4 contributions to the DY process as a function of the DY mass. We find 
that for the nominal LHC energy 14 TeV, the leading twist 2 contribution coincides with the all-twist result for masses 
down to about M = 6 GeV. Below that value resummation of the higher twists contributions is necessary. We have 
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a) b) c) d) 

FIG. 1: The leading and next-to-leading order diagrams for the Drell-Yan production. The diagrams c) and d) are enhanced 
in the small- a; limit due to a strongly rising gluon distribution. 



performed the analysis for both the transverse and longitudinal components and find that the twist 4 contributions 
have different signs in both cases (for ver low masses), however we do not observe the cancellations of the type found 
in the case of the DIS process fl7j . 



II. DRELL-YAN CROSS SECTION 



In the lowest approximation, the Drell-Yan lepton pair of mass M is produced form annihilation of two quarks 
of the same flavour / from the colliding hadrons: QfQf — > 7* — > see Figure QJi. In the collinear factorization 

approach, the leading order (LO) Drell-Yan cross section is given by 



d 2 a LO _ 47ra 2 m x\x 2 sr^ 2 
dM 2 dx F ~ 3N C M 4 xi+x 2 ^ 6 



}{q f ( Xl ,M 2 )q f (x 2 ,M 2 ) + q f (x ll M 2 )q f (x 2 ,M 2 )} , (1) 

/ 



where a em is the fine structure coupling constant, N c = 3 is the number of quark colors, qf/qf are quark/antiquark 
distributions in the colliding hadrons computed at the factorization scale /i 2 = M 2 and x\^ 2 are the light-cone momen- 
tum fractions of the quarks entering the scattering. In the LO approximation, the energy-momentum conservation at 
the photon vertex, {xyp + x 2 p) 2 — M 2 , leads to the following relation 

Xl x 2 = M 2 /s EE r , (2) 

where s = (p + p) 2 — 2p ■ p is the center-of-mass energy squared of the colliding hadrons. Introducing the Feynman 
variable of the lepton pair, xf = x\ — x 2 , one can easily find 



x i = \^ x f + 4r + x F ) , x 2 = \{yfxp + 4r - x F ) . (3) 

In the next-to-leading order (NLO) approximation, additional emission of a parton (quark or gluon) into the final 
state has to be taken into account. This is shown by the diagrams in Figure [Tp-[lJl- Because of the emission, the quark 
entering the photon vertex carries a fraction z < 1 of the original parton momentum. Thus, the energy-momentum 
conservation at the photon vertex, e.g. (x%p + zx 2 p) = M 2 , gives now 

xix 2 = t/z, (4) 

and the parton momentum fractions take the form 



xi = 5(\A4 + ^(rl z ) + x f) , x 2 = \ Wx 2 F + 4(t/z) - x F ) ■ (5) 

From the parton model conditions, xi. 2 < 1, we find that z > z m i n — r/(l — x F )- The NLO correction to the 
Drell-Yan cross section, proportional to the strong coupling constant a s , is given in the MS factorization scheme by 

m 

d 2 a NLO ijTa 2 m as (M 2 ) f 1 Xl X 2 ^ 2 



dM 2 dx F 3N C M 4 2tt 



f dz^j^Y^e}{q f { Xl ,M 2 )q f {x 2 ,M 2 )D q {z) 

j Zmin xi + x 2 t-r 1 1 1 



+ g(x u M 2 ) [q f (x 2 ,M 2 ) + q f (x 2 ,M 2 )] D g {z) + ( Xl f+ x 2 )} , (6) 
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FIG. 2: The DY cross section in the collinear approach with CTEQ6.6M parton distributions as a function of center-of-mass 
energy E — y/s at fixed xf = 0.15 and M = 8 GeV for pp and pp scattering. 



where x\ i2 are given by relations ([3]) and g is a gluon distribution. The coefficient functions D q g are of the form 

DM = cJ4(l + ^)(M^) + -2i^ln, + 5 (l-,)(|!-8)} , (7) 

DM = T fl |(z 2 + (l-z) 2 )lnii-^ + i+3z-Iz 2 | , (8) 

with Cf = 4/3, Tji = 1/2 and the standard (+) prescription to regularize soft gluon emission singularity at z = 1. 
Thus, the DY cross section in the collinear approximation up to 0(a s ) is the sum 

d 2 a c ° l d 2 a LO d 2 a NLO 

(9) 



dM 2 dx F dM 2 dx F dM 2 dx F 



The higher order (NNLO) corrections [SHlOj. proportional to a 2 , lead to much more complicated formulae for the 
coefficient functions and we will not provide them here. 



III. pp VERSUS pp SCATTERING 



We wish to compare the collinear factorization DY cross sections at the two presently operating colliders, the 
Femilab Tevatron which scatters pp and the CERN LHC which collides pp. Going from the proton into to the 
antiproton beam we have to interchange the quark distributions: u f-> u, d <-> d etc.. 

Such a replacement will not affect the singlet quark distribution in the second term in eq. ([7]) but it will modify the 
terms with the product of the quark distributions. For example, for two flavours, in the LO approximation 

(top® ~ e 2 l {u(xi)u(x 2 ) +u(xi)u(x 2 )} + e 2 d {d{x 1 )d(x 2 ) + d(xi)d(x 2 )} , (10) 

but in the pp case we find 

dftfp ~ e 2 u {u{ Xl )u{x 2 ) + u{ Xl )u{x 2 )} + e 2 d {d{ Xl )d{x 2 ) +d{ Xl )d{x 2 )} . (11) 
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FIG. 3: The enhanced diagrams from Fig. [T] in the proton rest frame. 



Both cross sections coincide if the momentum fraction x 2 is small. In such a case, to a good approximation, a sea 
quark is involved in the scattering since 

u(x 2 ) ^ u(x 2 ) , d(x 2 ) ^ d(x 2 ) ■ (12) 

This is illustrated in Fig. [2] where we show the DY cross section for pp and pp scattering as a function of the center- 
of-mass energy of colliding particles E at fixed xp = 0.15 and M = 8 GeV. For M <C E, x\ ss xp ~ 1 and 
x 2 w r /if <C 1, thus the two cross sections coincide. We use the NLO CTEQ6.6M parton distribution functions [2(| 
for this comparison. 



IV. SMALL x LIMIT 



The small x or high energy limit of the DY process means that dilepton mass is much smaller than the center-of-mass 
energy of colliding particles, M -C \fs- When x\ ~ 1 we have 

M 2 

x 2 = <1, (13) 

S Xl 

i.e. in the parton model, a fast quark (or antiquark) with the momentum fraction x\ annihilates with a slow antiquark 
(or quark) with the momentum fraction x 2 . In such a case, the diagrams in Figure [TJ;-[T}i, with a slow gluon, arc 
particularly enhanced due to the strongly rising gluon distribution in the small- X limit. Now, a difficult task of small- a; 
resummation arises (2lj . which touches the problem of unitarity corrections to the standard, linear QCD evolution 
equations. 

One can reformulate this problem in the rest frame of one of the protons, which acts as a target. In this frame, the 
diagrams in Figure [TJ; -QJl can be interpreted as the lowest order description of the process in which the fast quark 
scatters off a soft color field of the target with emission of a massive photon before or after the scattering, see Fig. |3] 
The photon subsequently decays into a pair of leptons. 

The cross section for radiation of a virtual photon from the fast quark of flavour /, which takes a fraction z of the 
radiating quark energy, is given by 12] 

Vt,l(iP^ 1* X ) = J d 2 rWl L (z : r : M 2 ,m f )a qq (x 2 ,zr) , (14) 

where T, L denotes virtual photon polarisation, transverse or longitudinal, respectively. Here r is the photon-quark 
transverse separation and w£ L are equal to [HI E3] 

W ( = ^ {[1 + (1 - z) 2 ] V 2 K 2 ( V r) + m 2 z^K^r)} , (15) 

W [ = ^M 2 (l-z)^W, (16) 

where -Ko,i are Bessel-McDonald functions, m/ is quark mass and rj 2 = (1 — z)M 2 + z 2 m 2 . 

The quantity a qq in eq. (|14[) is a dipole cross section known from DIS scattering at small Bjorken-x [22j . It was 
determined from fits to HERA data on the proton structure function F 2 at small x under different assumptions, e.g. 
assuming phenomenological form with a saturation scale Q 2 {x) — (x/xo)~ x [H, HI]: 

a qq (x,r) = a {l-eM-r 2 Q 2 s (x)/i)} . (17) 
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DY data from E772: E = 38.8 GeV 
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FIG. 4: The DY cross section from the collinear (solid lines) and dipole (dashed lines) approaches as a function of the dilepton 
mass M against the E772 collaboration data. The CTEQ6.6M parton distributions are used in the collinear case. 



Substituting this dipole cross section into eq. (fT4")l one can test predictions on the DY cross section in which parton 
saturation effects are taken into account. The final form of the DY cross section for the forward dilepton production 
is found after taking into account the incoming quark distribution in the proton 



dM 2 dx F ~ 6ttM 2 xi+x 2 ^f &f J Xl z 2 i qf 



(^,M 2 )+77 / (^,M 2 )]4 L (gp^ 7 *Y). (18) 



The expression in the squared brackets under the integral is proportional to the LO contribution of flavour / to the 
i<2 proton structure function, = e 2 x(qj + ?/)• Thus, the final formula reads 



a2„dy r i , 



dM 2 dx F 6ttM 2 



(19) 



A similar expression was found in by changing the variables r and z to p — zr and a = (1 — z)jz. The new 
variable p has the interpretation of a size of a quark pair while a is a quark/antiquark longitudinal momentum fraction 
with respect to the photon momentum. 



V. PREDICTIONS FOR THE LHC 



In Figure [4] we present a comparison of the results from the collinear factorization formula ([l} and the dipole 
formula (|19|) against the data from the Fermilab E772 collaboration (24j]. We use the next- to-leading order (NLQ) 
CTEQ6.6M parton distributions [13] for the collinear formula (solid lines) while the GBW parameterisation [25[ of 
the dipole cross section was used in the dipole formula. In eq. (flU)) . there are three massless light quarks and charm 
quark with mass m c = 1.4 GeV. For the energy y/s = 38.8 GeV and the indicated values of M and xf, the fraction 
of the slow parton momentum X2 ~ 0.01 — 0.1, which is slightly beyond the applicability of the dipole formula. We 
observe that the data are above the results from both the dipole and collinear factorization approaches. A similar 
result was found for the NLO MSTW08 parton distributions [26|. 
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FIG. 5: Left: the DY cross section from the collinear (solid lines) and dipole (dashed lines) approaches as a function of energy 
E — ^/s for fixed xf ~ 0.15 and dilepton mass M = 6, 8, 10 GeV (from top to bottom). Right: the same for M = 10 GeV and 
the indicated parameterizations of the parton distributions and dipole cross section. 



In Figure [5] (left) we present predictions for the DY cross section as a function of the center-of-mass energy E = ^fs 
at fixed xf = 0.15 and dilepton mass M = 6, 8, 10 GeV. For the collinear factorisation results we use the CTEQ6.6M 
parton distributions. In the figure on the right, we show the same results for M = 10 GeV in a more detailed way, 
using the linear scales. We additionally show the collinear factorization results for the MSTW08 parton distributions 



and the dipole approach results with two parameterizations of the dipole cross section: GBW [18( and GS 25[. In 
the latter parameterization the DGLAP evolution of the di pol e cross section for small dipole sizes is built in. We 
also analysed the Color Glass Condensate parametrization [27|, finding results very close to the GBW curves. At 
the LHC energy, the fraction X2 ~ 3 • 10 -6 and we are really in the small- x domain which has not been explored 
experimentally yet at the hard scale given by the invariant mass of the DY lepton pair. Thus, the presented results are 
only extrapolations. Nevertheless, we clearly see that saturation effects encoded in the dipole approach give results 
which are systematically below the collinear factorization predictions. At the LHC energy, in the most extreme case, 
the suppression of the DY cross section due to saturation effects can be as large as a factor of three. 
We are looking forward to the experimental verification of this result. 



VI. TWIST EXPANSION OF THE DIPOLE FORMULA 



We will analyze the dipole DY cross section (fT9| from the point of view of the twist expansion in positive powers 
of the ratio Q 2 (x2)/M 2 , where Q s is the x-dependent saturation scale. We utilize the Mellin transform, using the 
methods elaborated in [Tt], ■ The twist analysis for the Drell-Yan production is slightly more complicated than in 
the case of the structure functions. This is due to the fact that the integral in z cannot be performed analytically. 
This is because it involves the structure function ^(xi/z) which is given as a parametrization to the experimental 
data. In order to compute the twist expansion we will present two methods. In the first one we will expand the 
integrand around the z = 1 point and perform the integrals analytically. This will give us approximate, analytical 
expressions for higher twists. In the second method we leave the integrals in z and devise a method to extract the 
twist contributions numerically. This method gives exact results for the twist expansion. 

The standard definition of the Mellin transform of a function /(r 2 ) reads 



f°° dr 2 

m= / ^r(r 2 )-V(r 2 ), 
Jo r 



(20) 
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while the inverse transform is given by 



c-f-ioo 



where c is a real number which has to be taken in interval (a, b) such that the integral (|20|) is absolutely convergent 
for a < IZej < b. 

We start the twist analysis from the transverse part of the cross section, and we assume that quarks are massless. 
The part of the DY cross section which corresponds to the transverse polarization of the photon (fl9t reads explicitly 

d 2 *» Y a 2 em 1 fUz 



— F 2 (p-,M 2 } / dr 2 [l + (l-z) 2 ]M 2 (l-z)K 2 (MrVT^ l)a q g(x 2 ,zr). (22) 



dM 2 dx F 6tt 2 M 2 xi + x 2 
Using the above definitions of the Mellin transform, the dipole cross section can be written as 



= / -J-yr I —(r l2 y~<a qq (xy). (23) 



We will perform the twist analysis using the GBW parametrization of the dipole cross section (flTf which is given in 
a closed analytic form. The advantage of this is that we can perform the integral over the dipole sizes analytically. 
The Mellin transform of the GBW cross section reads 



r°° df 2 

<j G( 7 ) = a Q / — (f 2 )^(l - e"H = -a T(- 7 ) , (24) 
Jo r 

and therefore has single poles for all non-negative integer values of 7. 

Using representation (|2"3")l together with the explicit Mellin transform of the GBW cross section (I24[) , we can rewrite 
the above expression as 

d 2 a° Y _ a 2 em 1 f 1 dz 



/ ^F 2 (^,M 2 )[ dr 2 [l + (l-z) 2 ]M 2 (l-z)K 2 (MrVT~z 

Yx 2 J Xl z \ z J J 



dM 2 dx F 6tt 2 M 2 xi + . 

where the contour of integration over 7 is chosen so that the constant c satisfies < Rec < 1 and Imc = 0. We can 
now perform the integration over the dipole size r. To this aim we define 

S tM s r ^Ki^y . ^r( T )r( 1 + 7 ,r ( 2 + 7 ) 

Jo 21( 5 +7) 
The transverse part of the DY cross section reads therefore 



dM 2 dx F fo 2 M 2 Xl +x 2 J c _ loo 2m U > TU > \ AM 2 



dz „ fx-] . , \ ,~ so, ( z 2 \ (27) 



Clearly the poles in the 7 plane control the behavior in M 2 . We need to evaluate the integral over the longitudinal 
momentum fraction z. Let us introduce the following notation 

I T „{x u z, M 2 ) = - F 2 (^,M 2 ) [1 + (1 - zf]{z 2 y . (28) 
z V z J 

We first perform the integral over z analytically by expanding the above expression around z = 1. We assume that 
F 2 (xi/z, M 2 ) does not have any singularities at z = 1, which is corroborated by the expressions for the structure 
functions at these values of x% ~ 1. 

Expanding the integrand and integrating term by term one obtains, 



f l dzY,I [ T\{xi,z = ^M 2 ){l-zf f^-Y = J24 k \(x 1 ,z = l,M 2 )(l-x 1 ) 1 -'y+ k — - 
J ^ k>o V z s k>0 7 



*. • (29) 

7 + k 
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Inserting this expansion into eq. (|27|) one obtains a general expression which allows to extract the twists systematically 
in powers of (1 — x\) 

dAPd^ = 6tt 2 M 2 Xl + x 2 E J^ loo 2™ G W ffr W i-TTfc (1 "- Tl) U^l-xoJ ' (30) 



The twist expansion then corresponds to taking residues of the different poles in 7 that appear on the right-hand 
side of (|3"0|) . One can also rewrite it to expose a more general structure of this expansion 



-EE (l-^) fc (^) 7 'c fe , 7c (x 1 ,lng?/(M 2 (l-, 1 )),lnAf 2 ), (31) 

7c >l fc>l-7c 



where the hrst sum is performed over the poles in eq. ([3H)l . The coefficients Ck :lc depend on In Q 2 /(Af 2 (l — xx)) which 
reflects the fact that there can be multiple poles, and the InM 2 dependence comes from the possible dependence of 
F 2 on M 2 . Since F 2 is evaluated at large values of x\ we expect this dependence to be very mild. 



A. Twist 2 



To extract the twist 2 in powers of Q s (x 2 )/M we will consider the leading term in the expansion of the function 
It around z = 1. Taking this term, the integral over z reads 



dz I T ^ ) Jx 1 ,z = l,M 2 ) 



1 



J dzF 2 ( Xl ,M 2 ) 



F 2 (x 1> M 2 )(l-x 1 ) 1 -"<- . (32) 



7 

Therefore integral over z gives the single pole in the 7 plane. The function G( 7 ) = — T(— 7) also has a single pole, 
which together gives the double pole in 7 = 1. The higher order terms in expansion of ([28)) will contribute to terms 
which are suppressed by powers of (1 — x%) as is evident from eq. pip . We will discuss these later. 
Taking the first term in the expansion in 1 — z, (|32[) we obtain the approximate formula 



dM 2 dx F 6n 2 M 2 Xl +x 2 y Jc-too 27r * 1 -7 V 4M2 ( 1 - zi) 

The leading twist extraction amounts to taking the above formula and closing the contour to the right and taking the 
contribution from the pole at 7 — 1. The approximate expression for the leading twist is therefore 



(0) = a 2 em a F 2 (x u M 2 ) 2 Q 2 (x 2 ) 
T ' 2 6tt 2 M 2 xi + x 2 4M 2 



4 ,25 2 Q 2 Qr 2 ) 
-7B - 1 + 7^) " 7 In 7^727 



3 rv 2 y 3 4M 2 (l-xi) 



0(1- Xl ). (34) 



Note that there is an additional dependence on the mass M 2 in the function F 2 {x\, M 2 ). As 27 ~ 1 this dependence 
should be very mild and should not affect too much the twist expansion we are performing . We stress that this 
expression is an approximate one in the sense that there are subleading contributions coming from the other terms in 
expression (I28[) . They will not modify however the logarithmic term as they are regular as 7 = 1. More specifically 
the first subleading contribution to the leading twist comes from 



jf dz lW( Xl ,z = 1,M 2 ) (t-^) 7 = / dz(F 2 (x 1 ,M 2 )(l-2 1 ) + x 1 F^(x 1 ,M 2 )) (j^j 



7-1 



- [F 2 (x 1 ,M 2 )(l-2 1 ) + x 1 F^x 1 ,M 2 )](l-x 1 ) 2 ^j—. (35) 

This contribution is the leading one for the twist 4, but it also will affect the twist 2 as it multiplies the single pole 
in 7 = 1 present in G( 7 ). It gives the following correction to the leading expression for twist 2, eq. (|34l) : 

(1) _ a 2 em *o [F 2 {x 1 ,M 2 )-x 1 F'{x 1 ,M 2 )\ Q 2 (x 2 ) A 

At ' 2 ~~6^AP x7+x~2 {1 ' Xl) ^M 2 ~ X 3 ' {6b) 
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Twist-2 DY cross sections for x F = 0.15 
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FIG. 6: The twist-2 contributions, shown by the dashed lines, for the transverse and longitudinal DY cross section at the LHC 
energy y/s = 14 TeV. The solid lines show the all-twist results while the dotted lines correspond to the approximate twist-2 
results given by eqs. (fM)) and 



which is suppressed by one power of (1 — x\). 

The above-presented method allows to systematically compute the contributions to the given twist in powers of 
{l — x\). The effectiveness of this calculation depends however on the rate of the convergence of the series in (1 — x\). 
In our kinematics x\ is about 0.1 — 0.2, hence this convergence is rather slow. Below we will present another semi- 
analytical method which allows to compute the leading twist expression exactly. 

As the double pole contribution has been evaluated already in (|34[) . the remaining single pole contribution comes 
from the pole in G(j) multiplying the non-singular in 7 = 1 part of the integral over z. It suffices to take the 
regularized integral 

1,M 2 ) _ r l d ^ zF 2 (^,M 2 )(l + (1 - zf) - F 2 ( Xl ,M 2 ) 



xi 1 z J Xl 1 z 

which is evaluated at 7 — 1. The integrand in this expression is regular at z = 1 and thus the integral can be evaluated 
purely numerically. The correction to eq. (|34j) is therefore 

A ( fc >0) _ c^n 1 4 ( Ql{x 2 ) \ f 1 zF 2 (^,M*)(l + (l-z)*)-F 2 ( Xl ,M*) 



T < 2 " 6tt 2 M 2 xi + x 2 3 V 4Af 2 J J Xl 1-z 
and the complete, exact leading twist expression equals 

2 DY(t=2) 

d °T _ a (0) , A (fc>0) r „ q x 

A similar twist-2 analysis of the longitudinal part of the DY cross section is given in Appendix A. 

In Fig. [6] we show the comparison of the calculation of the all-twist formulae, eq. (|27|) for the transverse and 
eq. (|A2I) for the longitudinal parts (shown by the solid lines), with the approximate leading twist contribution. The 

approximate twist-2 contributions Ay 2 and A^, eqs. (|34|) and (|A6|) respectively, are shown by the dotted lines while 
the exact twist-2 expressions, eq. (|3"9")l and eq. (|A7I) . are depicted by the dashed lines. We see that the approximate 
twist-2 formulae are significantly above the exact result even for lager valaues of M. This is due to the large value 
of the subleading terms in the expansion of (1 — x\). Since x\ ~ xf = 0.15 the higher order terms can be still 
contributing. On the other hand, we observe that the exact twist-2 formulae are getting close to the all-twist result 
in the region M > 6. 
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FIG. 7: Convergence of the twist expansion of the transverse and longitudinal DY cross sections (left) and of their sum 
(right). The solid lines show the all-twist results while the dashed and dotted lines correspond to the twist-2 and twist-(2+4) 
contributions, respectively. 



B. Twist 4 

In this section we provide semi-analytical expressions for the twist 4 both for the longitudinal and transverse 
polarizations of the DY photon. We start with the transverse part and use eqs. ((27)) and (|35|) together, which will 
give part of the twist 4 contribution. Taking the residue at 7 = 2 we obtain 



(i) _ "L^o 1 ( ggfraA ± 
TA 6tt 2 M 2 x x + x 2 V 4A/ 2 J 15 



F 2 ( Xl ,M 2 ) ^-63 + 36 7jE + 18^(7/2) - 18 log 
+ x 1 F^(x 1 ,M 2 ) ( 17- 12 7B -6V(7/2) + 61o. 



V4M 2 (l-a;i) 

( ggCgg) 

^4M 2 (1 -xi) 



(40) 



The second part contributing to twist four comes from the contribution of the single pole in G( 7 ) multiplying the 
finite part in the integral over z. This can be evaluated in the similar manner as before by modifying the integral in 
order to perform analytical continuation to 7 = 2. To this aim we subtract in the integrand the first two terms of the 
expansion in (1 — z) 



ST(xi,M 2 ) 



1 ^ Z 3 F 2 (^,M 2 )(1 + (1 - zf) - F 2 (x u M 2 ) - (1 - z)[-3F 2 (x u M 2 ) + x^jx^M 2 )] 

(1-z) 2 ' 1 J 



which gives convergent integral by definition. Then we need to add the term coming from the lowest order term in 
expansion in z — 1 evaluated at 7 = 2 which is 

F 2 {x u M 2 ) 



1 — Xi 



The expression contributing to twist 4 is thus 

,,2 1 1 1 1 2 



a(2) ^m°0 1 ( y;(x 2 ) \ f _8 

T < 4 6tt 2 M 2 Xl + x 2 \ 4Af 2 J { 5> 



dF(xi,M ) 

1 — Xl 



(42) 
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where we evaluated the residue coming from the single pole at 7 = 2. The final expression for twist 4 for the transverse 
part is therefore 

= A (D , A (2) m 

dM 2 dx F ~ + T - 4 • [ 6) 

In Fig. [7] (left) we show convergence of the twist expansion for the transverse and longitudinal DY cross sections 
(left plot) and their sum (right plot). The solid lines show the all- twist results while the dashed and dotted lines 
correspond to the twist- 2 and twist- (2+4) contributions, respectively. We see that there is a poor convergence below 
M < 6 GeV. Both twist contributions are below the exact results for the transverse part of the cross section while 
for the longitudinal part they are above the exact result. This is the reason why the sum of the transverse and 
longitudinal parts of the twist contributions approaches the exact total cross section. This effect is shown in Fig. [7] 
(right). 



Summary 



In this paper we have investigated Drell-Yan production at forward rapidities and at high energies accessible at the 
LHC. We have used the dipole formulation suitable for forward rapidities and used the cross section which incorporates 
saturation effects. The comparison with the standard collinear formula shows the suppression of the production cross 
section for highest energies when using the dipole models with saturation effects. 

Using the dipole formulation with the GBW formula for the dipole model cross section, we have constructed the 
twist expansion for this process. Unlike the DIS case, where the twist expansion for the GBW model could be 
performed completely analytically, here we had to resort to the semi-analytical evaluation of the individual terms 
in the expansion. It was shown that the leading twist is a good approximation to the full result for masses of the 
Drell-Yan pair larger than about ~ 6 GeV. For lower masses the twist expansion quickly becomes divergent and full 
resummation is necessary. 
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Appendix A: Twist decomposition of the longitudinal part 



The leading twist-2 contribution for the longitudinal part of the cross section can be computed in a similar manner 
as for the transverse part. We use 



JO ^ i (2+T) 

After performing the integrals over the transverse coordinate and using the Mellin representation we obtain 



d °L __ "em°"0 



(Al) 



dM 2 dx F 6tt 2 M 2 xi + x 2 
Similarly to the transverse case, we can introduce the function 

I L (x u z,M 2 ) = - F 2 (^,M 2 ) (z 2 r . (A3) 
z V z ) 

The leading term coming from the expansion around z = 1 will thus give 

£dzl£\x 1 ,z = l,M 2 )(j^y = j\zF 2 (x 1 ,M 2 ) [j^y =F 2 (x 1 ,M 2 )(l-x 1 ) 2 -^^-^ . (A4) 
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Note that this expression is not singular at leading pole 7=1. Again the leading twist extraction in the longitudinal 
case therefore amounts to taking the formula 

,vr y„, !f 4^) (1 _ iil , r i c(i)ii(i) i fJ? Ly i 



dM 2 dx F 6ir 2 M 2 xi + x 2 y Jc-ioc 2m y " w/ 2 - 7 V4M 2 (1 - x x ) 

and closing the contour to the right, enclosing the pole at 7 = 1. The final result is 



{ o } ^ aj n <To2F 2 ( Xu M 2 ) Q 2 s {x 2 ) 2 _ 2 

AL ' 2 "6^Af2 z 1+X2 (1 Xl) 4Af 2 X 3 +U((1 Xl)) - {Ab) 



There are no logarithmic corrections because of the single leading pole coming only from the function G(j). 

Using analogous method as before we can also evaluate the exact expression for the longitudinal part of the cross 
section. The integral over z in (IA2|) is well defined for 7 = 1 and so to evaluate the exact twist 2 we can perform this 
integral directly numerically. The exact expression for the twist 2 contribution to the longitudinal part of the cross 
section reads 

d2 ° D S iT ; 2) = f&»*° 2 ggjg) x 2 [\ zzF2 (E1,m*) . (A7) 



dM 2 dx F 6ir 2 M 2 xi+x 2 4Af 2 3 



X\ ' Z 



Analogously we evaluate the twist-4 contribution to the longitudinal part. The first part of this contribution is 
obtained by taking the residue at 7 = 2 of expression (| A5|) 



(1) a 2 em a 2F 2 { Xl ,M 2 ) ( Q 2 s (x 2 ) V f 16 
L > 4 6^ 2 M 2 xi +x 2 \ AM 2 ) I 15 



(A8) 



The second part comes again from the finite part of integral in z which is 

r(o) 



1 I L , 7 =2(xi,z,M 2 )-I ( ">{xi,z=l,M 2 ) /-i z 3 F 2 (^,M 2 )-F 2 ( Xl ,M 



2 

X\ ^ " J x± 



r2\ 



dz " v "' 7= " v ' : i=/ rfz ^ ^ i. (A9) 

1 — z .L, 1 - z 



We thus obtain 



6tt 2 M 2 an + a; 2 \ 4Af 2 j \ 15/j,, 1 - 2 V ; 



The total contribution to the twist 4 for the longitudinal case is therefore 

2 W(t=4) 



do 



L _ A (l) , A (2 



dM 2 dxF 



= A^ 4 + A^. (All) 
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